Chemical pathways in ultracold reactions of SrF molecules 
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We present a theoretical investigation of the chemical reaction SrF + SrF — > products, focusing 
on reactions at ultralow temperatures. We find that bond swapping, SrF + SrF — > St2 + F2, 
is energetically forbidden at these temperatures. Rather, the only energetically allowed reaction 
is SrF + SrF — > SrF2 + Sr, and even then only singlet states of the SrF2 trimer can form. A 
calculation along a reduced reaction path demonstrates that this abstraction reaction is barrierless, 
and proceeds by one SrF molecule "handing off" a fluorine atom to the other molecule. 

PACS numbers: 34.20.Gj,34.50.Ez,34.50.Lf 



I. INTRODUCTION 

After many years of experimental effort, chemical re- 
action dynamics has now entered the cold and ultracold 
regime. Robust techniques such as buffer-gas cooling [l| 
and Stark deceleration |2j have pushed the energy res- 
olution of molecular beam techniques down to the mK 
level, resulting in new probes of chemical dynamics 

a- 

6]. Extending the low-energy limit even farther, coher- 
ent optical techniques have produced samples of alkali 
dimers in their absolute ground state, at temperatures 
in the 1-100 /^K range |7H10||. These ultracold molecules 
are so exquisitely sensitive to comparatively weak influ- 
ences, that chemistry can be studied and controlled by 
exploiting quantum statistics fill - electric fields [l2[ , and 
confinement in optical lattices |13l4l6| . 

With these new capabilities naturally come questions 
of what can be learned about chemical reaction dynam- 
ics under these novel circumstances. The manipulation 
of chemical reactions by external electric and magnetic 
fields relies heavily on the behavior of long-range physics, 
where the molecules exert, say, dipolar forces on one an- 
other, but are too far away from one another to react. In- 
deed, the theoretical analysis of the reaction 2KRb — > K2 
+ Rb 2 , observed and studied at JILA, took this point of 
view, by treating the actual reaction as a nearly perfectly 
absorbing "black box," which removed the molecules 
when they got close enough together, but without regard 
for what exactly happened to them [T3, EH ■ 

Vice versa, to understand more about chemical dynam- 
ics from these experiments would presumably require a 
scattering simulation on a complete four-body potential 
energy surface (PES) for the K-K-Rb-Rb system, which 
does not yet exist. In very recent work, however, Byrd et 
al. have described the main salient features of this sur- 
face [l9[ . First, along its main reaction coordinate, the 
reaction presents no energetic barrier to reaction (con- 
sistent with the high reaction rates observed at JILA). 
Second, its transition state, at the borderline between 



reactants and products, represents a T-geometry charac- 
teristic of an insertion reaction, wherein a K atom from 
one molecule inserts itself between the K and Rb of the 
other. This circumstance suggests that the reaction pro- 
ceeds by a complicated four-body dance of the involved 
atoms, which may be revealed at ultralow temperatures 
by characterizing the resonant states of the complex. Ob- 
serving these transition-state resonances and interpreting 
their influence in chemical reactions is a longstanding 
goal of physical chemistry, and one in which the high 
energy resolution of ultracold molecules may be a great 
help. 

Set against this backdrop, it would also be useful to ex- 
plore other molecules reacting at ultralow temperatures, 
to gain additional insight into what can be learned. To 
this end, the SrF molecule is an appealing candidate, and 
the one with which we deal in this article. SrF is a prime 
example of a class of molecules that have been identified 
as amenable to direct laser cooling from a beam [20|, and 
in fact laser cooling has been recently demonstrated [2l[ . 
It is quite polar (dipole moment » 1,4 Debye), so that 
electric field manipulation is a possibility. Moreover, it 
has an open-shell 2 £ ground state that gives it a magnetic 
moment as well. This circumstance opens opportunities 
such as trapping the molecules magnetically, while ma- 
nipulating their interactions electrically (22j . 

In this article we explore the possibility and mecha- 
nisms for chemical reactions of SrF molecules at ultralow 
temperatures. We make several observations. First, the 
exchange reaction 



2 SrF -> Sr 2 + F 2 + AE e 



(1) 



is energetically disallowed, as might be expected for a re- 
action that turns two ionic bonds into two covalent ones. 
Therefore, the reaction, if it happens at all, must pro- 
ceed by an abstraction reaction in which an atom jumps 
from one molecule to the other. This is exactly the op- 
posite situation from what occurs in KRb, where the ex- 
change is the only possible reaction, and then only just 
barely 0, H Ej 

Second, we find that the Sr-abstraction reaction 
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2SrF 



Sr 2 F + F + A^ rlmel 



(2) 



2 



cannot occur at low temperature, whereas the F- 
abstraction 

2SrF ->■ SrF 2 + Sr + AE ttimel (3) 

can occur, as it produces the deeply bound SrF 2 trimer. 
Moreover, this reaction, which is barrierless, can only 
occur in the singlet channel, i.e., only if the reaction 
takes place on the PES with total electronic spin 5 = 0. 
Therefore, chemical reactions are expected to occur at 
quite high rates for unpolarized SrF molecules, while 
they should be strongly suppressed for spin-polarized 
SrF, which scatter primarily on the triplet surface. Spin- 
rotation couplings will ensure that this suppression is not 
complete |25|]. Nevertheless, one must be mindful of any 
opportunity to suppress inelastic collisions, as they can 
easily destabilize the gas and derail attempts to exploit 
the molecules for many-body physics applications. 

Third, we describe the basic physics of the abstraction 
process by looking at a restricted version of the four- 
body PES. We find that the F end of one SrF molecule 
approaches the Sr end of the other, as dictated by the 
dipole-dipole interaction between molecules. Next, the F 
atom is handed off from one molecule to the other, and 
the free Sr goes off by itself. By this hand-off mechanism 
it is possible that the reaction proceeds without forming a 
resonant complex, and that its interpretation from ultra- 
cold collision data may be more straightforward than the 
complete re-arrangement necessary in reactions of KRb. 
In any event, complementary dynamics should give us 
complementary insights into these quite different reac- 
tions. 



II. AB INITIO CALCULATIONS 

In this section we will study the relevant molecular 
species: the dimers SrF, Sr 2 and F 2 ; and the trimers SrF 2 
and Sr 2 F. All calculations of these species are performed 
using the molpro suite of ab initio electronic structure 
codes [26| . We use the relativistic effective core potential 
and associated basis set of the Stuttgart group for Sr [27j 
(ECP28MDF) and the augmented, correlation consistent 
valence triple (quadruple) zeta basis set (AVTZ (AVQZ)) 
of Dunning for F (28[. For the diatomic properties, the 
AVQZ basis is used. 

We perform the calculations by first computing a spin- 
restricted Hartree-Fock (RHF) wave function as a start- 
ing guess for a coupled-cluster singles, doubles, and non- 
iterative triples excitations calculation (CCSD(T)) [29(. 
The minimum energy configuration is obtained using the 
method of steepest descents. Properties such as dipole 
moment, polarizability, and quadrupole moment are cal- 
culated using the finite field approach. To compare the 
energies to the free atom limit, we use basis set superpo- 
sition error (BSSE) corrections as given by the method of 
Boys and Bernardi [30] ■ We calculate the vibrational fre- 
quencies by making the symmetry-independent displace- 



ments of the atoms and calculating the Hessian within 
the molpro suite of routines. 



A. Diatomic species 

In the collision of two SrF molecules there are three 
diatomic species to consider; the reactant SrF and the 
possible products Sr 2 and F 2 . In the following we will 
investigate the relevant diatomic properties such as equi- 
librium bond length, well depth, and vibration constant. 
But first, an exercise in bond strengths can immediately 
inform us as to whether the reactions in Eq. ([!} can oc- 
cur. On the left hand side of the reaction we have two 
ionic bonds. On the right hand side we have one van 
der Waals bond and one covalent bond; the former is 
very weak and the latter usually less bound than an ionic 
bond. An ionic bond can be thought of as a strengthen- 
ing of the single covalent bond because of a transferring 
of charge from one atom to the other. From these sim- 
ple considerations, we conclude that the reaction in ([T]) 
is energetically unfavorable at ultracold temperatures. 

Sr 2 is a van der Waals molecule. The 1 So atomic struc- 
ture means that the the outer s-electrons are already 
paired up and therefore play a limited role in the bond- 
ing of two Sr atoms. Previous work has shown that this 
molecule is not deep ly bound, with a well depth of D e — 
1081.8 cm -1 [UHl]. Because of the large binding length 
and relatively larger mass of Sr atoms, the vibrational 
constant is fairly small, ui e — 40.3 cm -1 [32J. Therefore, 
the energy required to dissociate this molecule from the 
zeroth vibrational level is given by Dq = 1061.6 cm . 

F 2 is a covalently bonded species, with a ground state 
of symmetry. F 2 has an appreciably large bind- 
ing energy of Dq =12950 cm -1 33] and much smaller 
bond length compared to Sr 2 . The unpaired p-electrons 
in the F atoms pair up and form a fairly deep well. 
Due to the lighter mass and tighter confining poten- 
tial, the F 2 molecule has a comparatively large vibra- 
tional constant, 917 cm -1 (34[. Therefore, the well depth, 
D e = 13410 cm -1 when combining the work of [33| and 

m 

SrF is a highly polar molecule. The willingness of F to 
take an extra electron, and of Sr to give one up, would 
lead one to conclude that at the minimum of the well the 
molecule is well described by a Sr + F~ configuration. We 
therefore compared the optimized geometry calculated at 
the RHF-CCSD(T)+BSSE calculation in two ways. The 
first was to compare the energy at the bottom of the po- 
tential energy surface to that of free Sr and F atoms. The 
other was to compare to free Sr + and F _ ions and then 
calculate the ionization potential of Sr and the electron 
affinity of F. These calculations agreed with each other 
to within several cm -1 . The well depth D e we report 
here refers to the energy required to separate SrF into its 
neutral partners, Sr and F. 

In Table Q] we present the results of our RHF- 
CCSD(T)+BSSE calculations. All bond lengths (R e ) are 
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TABLE I: The molecular properties of Sr2, F2, and SrF. En- 
ergies are in cm -1 , bond lengths in A, and moments in atomic 
units. Experimental values are given where known. 



Molecule 


D e 


Re 


We 


dm 


a z 


6» zz 


SrF (this work) 


44200 


2.084 


499 


1.38 


126 


8.95 


Expt. 


45290(560)" 


2.075 6 


502.4 6 


1.36 c 






Theory 


45000 d 


2.085 


507 








Sr2 (this work) 


820 


4.773 


36.2 








Expt. E 


1081.8 


4.672 


40.3 








F2 (this work) 


12880 


1.410 


927 








Expt. 


13410 / 


1.411 9 


917 9 









"Ref. 
'Ref. 
c Ref. 
d Ref. 
e Ref. 
fRef. 
9Rof. 



35]. 
36]. 
37]. 
38]. 
32 I . 

32 I, ZEKE/ion-pair imaging Do adjusted by co e of [34 

a. 



e 



FIG. 1: (Color Online) Diagram of the geometries considered 
in the calculation of the the molecules Sr2F and SrF2. The 
coordinates ri, ri, and 9 were varied and energy optimized 
to find the minimum energy configuration in both the A' and 
A" symmetry groups. 



in A, well depths (D e ) and vibration constants (ui e ) are 
in cm , and associated dipole and quadrupole moments 
(d m and 9 ZZ ) are in atomic units (e ao and e Oq, re- 
spectively). The corresponding experimental values are 
from the references given, and uncertainties are printed 
for those where the value was reported. The theoret- 
ical calculation of Langhoff et ai, is done at the con- 
figuration interaction with single and double excitations 
level of theory [38]. Our current method yields a well 
depth smaller than the experimentally obtained value by 
2(j I35j . The bond length is in good agreement with that 
of [23]. From these results we immediately see that the 
reaction in Eq. ([T]) is energetically forbidden, by more 
than 70,000 cm -1 , and is of no concern in an ultracold 
gas. 



B. Triatomic species 

There are two different triatomic species to consider 
as given in Eqs. (J2J) and ([3]). At first glance, we might 
think both are quite easily allowed. From a simple bond 
strength argument, we can see that one is favored over 
the other. Particularly, that the triatom in Eq. @ is less 
bound than the triatom in Eq. ^ . The triatom in ([3]) has 
two ionic bonds, while the triatom in ^ can have only 
one. Therefore, from the fact that two ionic bonds should 
be more deeply bound than one, we expect the triatom 
in (J3j) to be more deeply bound than in {2j). In addition, 
it is expected that the reaction in @ is energetically 
unfavorable because there are stronger bonds on the LHS 
than on the RHS. However, the reaction in §3§ cannot 
be determined from these arguments. Thus, it is to be 
calculated using the same ah initio methods as in the 
previous section. For the SrF2 singlet state, we use the 
AVQZ basis for F as we did for SrF. However, for the 
remaining states we only used the AVTZ basis. This is 
because these states are so energetically forbidden (as 



TABLE II: The molecular properties of SrF 2 and Sr 2 F. En- 
ergies are in cm -1 , bond lengths in A and angles in degrees. 
Previous theoretical values are given where applicable. 



Molecule 


Symm. 


D e 


ri 


ri 


9 


SrF 2 
Theory" 
Expt. 6 


X A: 


91870 
89300 to 91900 


2.13 
2.16 


2.13 
2.16 


135.7 
138.8 


Sr 2 F 


' 2 B 2 


52300 


2.26 


2.26 


115.8 


Sr 2 F 




43900 


2.25 


2.25 


101.7 


SrF 2 


J A" 


43400 


2.10 


3.72 


70.4 



"Ref. & 

'Ref. [401 ] . Extracted from enthalpy of formation. 



we will see) that the less time-consuming calculation is 
adequate for our purposes. 

In Fig. 1, we give the three coordinates that are used in 
the optimization of the geometry. Using these three co- 
ordinates (r%, r2, and ff), without constraining then, en- 
sures that the electronic wave function belongs to the C s 
symmetry group. Thus, there are electronic states with 
even (A') and odd (A") reflection symmetry through the 
plane containing the three atoms. When the two bond 
lengths are equal, the symmetry group is instead ex- 
panded to C2V This added symmetry describes whether 
the molecule possesses even or odd reflection about the 
line which bisects the homonuclear bond. Where the 
ground state has equal bond lengths, we have indicated 
this with the appropriate notation; A\ (B2) have even 
reflection symmetry about the mid-point between the 
homonuclear bond and even (odd) reflection symmetry 
through the plane containing the three atoms. Similarly, 
A2 (Bi) has odd reflection about the midpoint of the 
homonuclear bond with even (odd) reflection through the 
plane containing the three atoms. 

In Table [H] we present the results of the RHF- 
CCSD(T)+BSSE calculation. As is evident, the most 
deeply bound trimer is given by the 1 K\ symmetry for 
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SrF2. It is also the only one for which previous the- 
oretical work is available for comparison. In Ref. [39j . 
Kaupp et al. present a theoretical study of many alkaline- 
earth metal atoms with two halogen atoms. We note that 
alkaline-earth di-halidcs have historically been of signifi- 
cant interest in chemistry, as some have bent geometries 
(hence possess permanent dipole moments in their body 
frame) and some are linear (with no such dipole moment). 
A bend is not expected from classical models of charged 
particles interacting, and therefore correctly predicting 
the bend is a measure of the quality of the wave func- 
tion. This is because the bend is a result of the presence 
of (^-electron orbitals on the ground state. Our calcula- 
tion agrees well with that of [39| in the description of the 
bond length and bending angle. For the first time, we 
present the energy required to pull apart the molecule 
into its atomic constituents - atomization. This energy 
is greater than that of two SrF molecules and therefore, 
in the absence of a barrier, should react chemically at 
ultracold temperatures. 

Vice versa, the triplet SrF2 molecule (last line of Ta- 
ble II) is far less bound, and indeed has binding energy 
similar to that of a single SrF molecule. A simple way 
to understand this is to think of this as a triplet covalent 
bond in F2 with a Sr atom bound to one or the other F 
atoms ionically. The triplet bond in F2 is much shallower 
than the singlet covalent bond in the ground state. Then 
the Sr comes along and binds to one of the F atoms. A 
simple charge population analysis shows that only the 
F atom closer to Sr captures charge. Importantly, the 
triplet trimer is energetically disallowed at ultracold (and 
even room) temperature studies of SrF collisions. 

The other potential trimer, S^F, is reported in the sec- 
ond and third lines of Table ITll This molecule, while more 
deeply bound than the triplet state of SrF2, is not deep 
enough to be chemically reactive, as expected from the 
simple bond arguments. Thus, this reaction will not pro- 
ceed at ultracold temperatures. Notice that this molecule 
has a binding energy only slightly more than that of a 
single SrF diatom. The doublet nature of this molecule 
is due to the one unpaired electron that can be viewed 
as coming from the F atom plus a Sr2 bond, or from the 
doublet nature of SrF plus a free Sr atom. 

We therefore conclude that only the singlet electronic 
state of SrF 2 will form in the collision of two SrF 
molecules, at ultralow collision energy. We find that the 
energy released is 3470 cm . To be thorough, we must 
take into account the zero-point vibrational energy of the 
systems in consideration. In Table U we gave the cal- 
culated and experimental values of the vibrational con- 
stants of SrF. We have calculated the vibrational con- 
stants of SrF2 in the singlet electronic state: u; sym = 
472.9, Wasym = 478.7, and Wbend = 79.3 cm— 1, where 
these represent the symmetric (or breathing) mode, the 
anti-symmetric mode, and the bending mode of the tri- 
atomic molecule. This yields a zero-point energy of 
515.5 cm -1 . Therefore, all the exoergic reactions are 
reduced by sa 20 cm -1 . It is evident that the added vi- 



brational modes (three vs one in the diatom SrF) are not 
enough to overcome the energy released in the chemical 
reaction. In fact, there are roughly 150 vibrational states 
accessible in the reactants, given using the calculated ex- 
oergicity of A tr imei = 3450 cm -1 . 

To our knowledge, this quantity has been measured 
only once, in a gas cell which maintained chemical equi- 
librium between reactants and products at ~ 1500 K, 
using mass spectroscopy to measure their relative abun- 
dance 41|. Knowing the temperature, this measurement 
provided information on the enthalpy of formation A.fH, 
which corresponds to the energy released, Ai5 tr i mer in 
Eq. ([3]). This measurement agrees that the reaction is 
exothermic, but by a more modest value of AfH = 2.1 
kcal/mol = 740 cm^ 1 , than we have calculated. Given 
that the experimental uncertainty was comparable to 700 
cm -1 and our computational uncertainty is perhaps 2000 
cm , the results are not too seriously in disagreement. 
Ultimately, measurements in laser-cooled SrF samples 
should sort out this issue in detail. 



III. REACTION PATH BASICS 

To compute the full PES of the four-body system in- 
volved in the reaction is of course a complicated affair. 
Here we instead compute selected slices of this PES. One 
is a reaction-path coordinate version that will verify that 
the reaction is barrierless. The second is a more detailed 
examination of the "hand-off" mechanism that drives this 
abstraction reaction. We will focus exclusively on the 
singlet PES, as it is the only one leading to reactions in 
an ultralow temperature SrF gas. All calculations are 
performed with the ECP28MDF basis set and ECP of 
the Stuttgart group for Sr along with the AVTZ basis of 
Dunning for F. We use the RHF-CCSD(T)+BSSE level 
of theory. 

A basic idea in constructing the reaction path surface 
is that the relatively heavy Sr atoms move comparatively 
slowly, and therefore the Sr-Sr distance can be regarded 
(approximately) as an adiabatic coordinate. Fixing this 
distance, denoted i?sr-Sr, we optimize the coordinates 
T\, 6i, and 82 of the F atoms, as defined in Figure 
2, so as to minimize the energy. For simplicity, we con- 
strain all four atoms to lie in a plane. While this limits 
the possibility of the molecules changing the dihedral an- 
gle, it does not prevent all insertion type configurations 
from being explored. As we will see, the optimizations 
of the other coordinates suggest a preference for linear 
orientations upon approach. 

To track the progress of the reaction, we define the 
approximate reaction coordinate A via 



r-2 



(4) 



In the limit that A is negative and large, it denotes the 
distance between the two SrF diatoms, whereas when A 
is positive and large, it refers to the distance between the 
SrF2 trimer and the free Sr atom. 
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FIG. 2: (Color Online) The top portion of the figure describes 
the geometry of the situation, n, r2, 0i, and 02 are optimized 
so as to produce the minimal energy for a fixed R$ x s t . The 
LHS describes the approach of two dipole objects. On the 
RHS, the approach is more van der Waals like where the polar 
SrF2 polarizes the Sr atom. The middle is a sharp point 
that is an artifact of the choice of A = r-2 — r' . The larger, 
blue (smaller, red) filled disks represent Sr (F) atoms. The 
approach to the transition state near A = A is very nearly 
linear up to the point of F capture. Then 6\ and 62 rotate 
to form a complex which contains a slightly distorted SrF2 
configuration. The lowest panel shows the optimized angles 
9i (solid, black line), 82 (large dashed, blue line), and 9 (small 
dashed, red line). The values remain fairly constant for |A| > 
1 A, and change to form the transition state in the region 
|A| < 1 A 



The optimized energy is plotted as a function of A 
in the middle panel of Fig. 2, adjusted so that the zero 
of energy refers to two free, separated SrF molecules. 
The PES along this path shows very little structure: it 
corresponds to almost pure dipolar attraction between 
reactants for A < 0, and almost pure van der Waals 
attraction between the products for A > 0. In particular, 
this cut through the PES shows clearly that there is no 
barrier to reaction, at least along the minimum energy 
path. 

This choice of reaction coordinate is useful, but also 
leads to a structural discontinuity in the surface of Fig. 2. 
This discontinuity is made clear in the last panel of Fig. 2 



which depicts the variation in the angles 9\, 62, and 
9 (see top panel of Fig. 2). This is an artifact of the 
choice in approximate reaction coordinate. Physically, 
this discontinuity arises because i?s r -Sr is not strictly an 
adiabatic coordinate. At some point, the intermediate 
F atom glides from its local minimum (attached to the 
right-hand Sr atom) to its global minimum (attached to 
the SrF on the left). When this occurs, a single value of 
^Sr-Sr corresponds to two distinct values of the reaction 
coordinate A where the adiabatic energy is the same. 
This mechanism will be clarified by looking at the PES 
from a different perspective, below. 

For A < 0, the system is described by a linear geometry 
(61 = 62 = 0°). As the diatoms approach one another, 
they reach a location in A where the discontinuity occurs. 
From the A > side, the values of 9i, 92, and 9 are fairly 
constant. As A is decreased toward zero by decreasing 
Rs r _sr, the 4-atom system reaches a minimum near A = 
0. As the Rsr- Sr is decreased further, 9\ approaches 0° 
and the system approaches the geometrical discontinuity 
from the other side. This leads to a fairly smooth reaction 
path in the coordinate A, but hides the change in the 
geometrical configuration of the 4-atom system. 

Using the information contained in the values of 9±, 
92, and 9 we are led to study the F capture process by 
analyzing the linear configuration. Even though the fi- 
nal trimer's ground state is bent, the energy of the linear 
configuration is only ~500 cm -1 above the bent configu- 
ration of SrF2. For this reason, we can make a reasonable 
qualitative description of the reaction by constraining all 
atoms to lie on a line, as in the upper panel of Fig. 3. 

We again treat the Sr-Sr distance as an adiabatic co- 
ordinate. Then, for fixed Sr-Sr distances i?sr-Sr, we vary 
the location of the middle F atom, describing it with 
coordinate r' as depicted at the top of Fig. 3. The left- 
most F atom is fixed to the diatom bond length of SrF, 
r\. Because the bond length in SrF and SrF 2 are so 
similar, this is a reasonable assumption. Using the RHF- 
CCSD(T)+BSSE method we calculate the potential ex- 
perienced by the middle F atom. The results are pre- 
sented in Fig. 3 for various values of i?s r -Sr- The BSSE 
compares the energy to the atomization limit. In the fig- 
ure, we offset the energy by the binding energy of two 
SrF molecules. Thus, zero energy is the energy of two 
SrF diatoms. 

In the first plot of Fig. 3 the two Sr atoms are far 
from one another. The potential for the F atom has 
two minima, representing the fluorine attached to the 
original Sr atom (right hand well) or else attached to the 
other SrF to form the product SrF 2 (left hand well) . The 
latter well being deeper signifies that the products SrF 2 + 
Sr are energetically favored and the reaction is exoergic, 
even within a strictly linear geometry. Between these 
minima stands a high barrier, which naturally prevents 
the F atoms from jumping from reactants to products 
when the Sr atoms are this far apart. 

In subsequent panels of Fig. 3, the two Sr atoms ap- 
proach one another. In each case the barrier lowers until 
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FIG. 3: (Color Online) The top portion gives the geometry of 
the system. The bond length ri is fixed while r2 is varied, r' 
gives the distance to the middle F to the other Sr atom and is 
called the hand-off coordinate since it describes the handing 
off of a F to one Sr from the other Sr. Each subsequent 
plot is for a fixed value of i?sr-Sr given in the upper RHS. 
As the value of i?sr-Sr is reduced, the barrier to forming the 
system F-Sr-F + Sr is diminished until it no longer stops the 
transfer of F. The calculations were performed at the RHF- 
CCSD(T)+BSSE with the Sr ECP28MDF and F AVTZ bases 
respectively. 



it eventually disappears altogether, and the F atom rests 
at the bottom of a single minimum in the last plot of 
Fig. 3. This configuration defines the four-body transi- 
tion state with A ss A. At this point the Sr atom can 
recede, and to the extent that its motion really is adia- 
batic, the F atom will tend to remain in the lower well, 
thus finding itself a part of the SrF2 final product. We say 
that this F-abstraction reaction occurs by a "hand-off" 
mechanism, whereby the right-hand Sr atom approaches, 
gently hands off the F, and then goes away. As compared 
to the KRb-KRb surface [l^| , this reaction is less likely to 
partake in a complicated dance of the four atoms, becom- 
ing thoroughly enmeshed in a four-body transition-state 
complex. 



IV. CONCLUSIONS 



In summary, we have established that ground-state 
SrF molecules will indeed by chemically reactive, even 
at ultralow temperatures that will be achieved via laser 
cooling. The only possible outcome of such a reactive 
collision would be the singlet trimer SrF2, plus a free 
Sr atom. Because of the need to produce a singlet final 
state, chemical reactivity should be strongly suppressed 
in a spin-polarized sample. This species therefore looks 
like a promising candidate for ultracold molecular stud- 
ies. On the one hand, SrF may serve as a useful platform 
for probing and understanding abstraction reactions in 
unprecedented detail. On the other hand, if producing 
a spin-polarized sample suppresses reactions sufficiently, 
there is also hope that the molecules may live long enough 
to perform interesting experiments on dipolar degenerate 
quantum gases. To estimate if this is so, future work will 
need to study in more detail both the triplet surface, and 
its detailed coupling to the reactive singlet surface. 
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